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Abstract 






The SNO result on charged current deuteron disintegration, the SuperKamiokande 1258- 
day data on electron scattering, and other solar neutrino results are used to revisit the model of 
T— I ■ neutrino oscillations driven by a violation of the equivalence principle. We use a x^ minimiza- 
^ . tion technique to examine oscillation between the fe and another active neutrino, both massless, 
pg , and find that within the Standard Solar Model the fit to the SNO and SuperKamiokande spec- 

^^ ' tra are moderately good while a very good fit is obtained when the absolute normalizations of 

C — ' s 

f— N ' the B and hep neutrino fiuxes are allowed to vary. The best fit prefers large, but not maxi- 
mal, mixing, essentially no hep neutrinos, and a 40% reduction in the ^B neutrino flux. The 
^^ • fit to the total rates from the different experiments is not encouraging but when the rates and 
"^ , spectra are considerd together the situation is much improved. We remark on the expectations 

I ' of the VEP model for the neutral current measurements at SNO. 

Oh! 
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1 Introduction 

The Sudbury Neutrino Observatory (SNO) collaboration has presented results from 241 
days running of their solar neutrino experiment |l|], a reminder that it is still early days 
yet in this fruitful area of astro-particle physics. From the data on charged current (CC) 
deuteron disintegration {ug + d —* p + p + e~) they observe a Ue flux which is 3.3a less 
than the neutrino flux measured by the Super-Kamiokande collaboration (SK) via u — e^ 
scattering (ES) 0Q. This further strengthens the oscillation interpretation of the solar 

*E-niail address: amitava@cubnib.ernet.in 
'E-mail address: arun@cubnib.ernet.in 

^SNO also has results from ES events ^ which agree with the SK results [|| but is of poorer statistics 
currently. 



neutrino problem and, at least in a two neutrino scenario, disfavours the mixing of the 
z/e with a sterile species. 

Violation of the Equivalence Principle (VEP) has been examined in the literature as a 
possible solution to the solar neutrino problem 0, ^, |^, P, |^ , one of the alternatives to the 
mass-mixing vacuum oscillation and the matter induced Mikheyev-Smirnov-Wolfenstein 
(MSW) approaches [§. In this note we revisit the VEP solution in a two- flavor picture 
with the z/e oscillating to another active species both of which are assumed to be massless. 

Neutrino oscillation in VEP can occur if (a) the weak equivalence principle is not 
satisfied, i.e., the coupling of neutrinos to the gravitational field is nonuniversal and 
(b) the flavor eigenstates are not identical to the states that couple to gravity 0. It 
does not require neutrinos to carry a non-zero mass. The weak equivalence principle 
demands the coupling of particles to the ambient gravitational potential (f) to be of the 
form f(j)E, where E is the particle energy, and / a universal coupling constant. If the 
latter varies from one neutrino species to another then that would constitute a violation 
of the equivalence principle. If /i 7^ /2 in a two-neutrino framework, then these states 
define a basis in the two-dimensional space which, in general, could be different from 
the flavor basis. 

Pe = ^1 cos 6 + 1/2 sin 6] v^ = —ui sin 6 + 1/2 cos 6, (1) 

where x = fi, t. The effect due to a small splitting A/ will manifest itself in the form 
of flavor oscillations, the wavelength going to infinity as A/ tends to zero. The essential 
difference between this approach and the mass-mixing solution appears in the energy 
dependence of the survival probability. For a two-neutrino picture, the general expression 
for the survival probability for an initial z/g after propagation through a distance L is 
given by: 

Pee{E,, L) = l- sin^ 29 sin^ (^^ , (2) 

where, for the VEP induced case, the oscillation length A is: 

Ey is the neutrino energy. Thus, unlike the mass- mixing case where X oc E,y, for VEP 
A oc l/E^. Due to the different energy dependences of the survival probability in the 
mass-mixing and the VEP alternatives, their predictions can be quite different for solar 
neutrinos of different energies. We follow the prevalent practice of choosing the gravi- 
tational potential, 0, to be a constant over the neutrino path. This is the case if the 



potential due to the Great Attractor |T^ dominates over that due to the sun and other 
heavenly bodies in our neighborhood. In such an event, writing AF = (pAf /2, the 
expression for the oscillation wavelength, Eq. (^, becomes 

, 6.20 X lO-^^m /lMeV\ 



An alternate physics scenario which also leads to neutrino oscillations with A oc 1/ Ey 
is a picture of violation of special relativity (VSR) |Tl| . If special relativity is violated, 
the maximum attainable speed of a particle in vacuo need not universally be the speed 
of light c. In particular, if the maximum possible velocities of the two types of neutrinos 
be f 1 and V2 and these velocity eigenstate neutrinos be related to the z/g and u^ through 
a mixing angle 6 (see Eq.(|l])) then in this case the expression for A is: 

where Av is the velocity difference for the neutrinos ui and 1^2- Comparing Eqs. (^ 
and (0) one finds that the energy dependence of the oscillation length is identical in the 
two cases0 and the role of Av in the VSR case is the same as that of 4>Af in the VEP 
formalism. Here, we use the terminology of the VEP mechanism but the results can be 
taken over mutatis mutandis to the VSR situation. 

In this work we use the BBP2000 calculation [|13| of the solar neutrino flux as the SSM 



reference. In addition to this, we explore the possibility of the absolute normalizations 
of the solar ^B- and /lep-neutrino spectraQ being different from their BBP2000 SSM 
predictions. If Xb and X^^p denote the factors by which the absolute normalizations are 
multiplied, we use the data to find the best-fit values for these. We find that all the fits 
are improved in a noteworthy manner when Xb is permitted to be less than unity. 

2 Update of VEP for the SK ES results 

In contrast to other recent work on VEP |^, ^] based on the 825-day SK spectral data, 
in a previous publication 0, the SK ES results from 1117 days of running |1|] were 



considered. It was found that the predictions were by- and- large robust. One trend, 
which was remarked upon, was that the softening of the highest energy SK ES data had 
resulted in a reduction of Xh^p- In fact, it was noted that the best-fit to the day-night 
spectrum of the 1117-day SK data preferred a negligibly small value of Xhep- 

For a further update, we have performed a similar analysis with the SK 1258-day 
data. The SK results are presented in the form of number of events (normalized to the 
SSM expectation) in 19 electron recoil energy bins of width 0.5 MeV in the range 5 MeV 
to 14 MeV and a 19th bin which covers the events in the range 14 to 20 MeV |@]. 

The definition of x^? the error correlations, etc. are chosen in the same manner 
as in 0]. Suffice it to say that we have included the statistical error, the uncorrelated 



systematic errors, and the energy-bin-correlated experimental errors ]T5|] as well as those 



from the calculation of the shape of the expected spectrum [0. When we allow the 



^It has been shown that inclusion of CPT-violating interactions in addition to Violation of Special 
Relativity can lead to more general energy dependences involving l/E^, E^, and constant terms |12| 
■^The SNO and SK experiments are sensitive only to the ®B and hep neutrinos. 



normalizations of the ^B and hep fluxes to vary, we do not include their astrophysical 
uncertainties separately. The results are presented in Table 1. 



Fitted 
Data 


Case 


sin^ 29 


AF 

(10-24) 


Xb 


-^hep 


xLn/d.O.f. 


g.o.f 

(%) 


SK1258 ES 

Day-Night averaged 

Spectrum 


la 


0.72 


0.22 


1.0 
(fixed) 


1.0 

(fixed) 


29.10/17 


3.37 


lb 


1.0 


0.99 


0.77 


2.4 X 10-^ 


11.38/15 


72.52 


SK1258 ES 

separate Day & Night 

Spectra 


Ic 


0.72 


0.21 


1.0 

(fixed) 


1.0 

(fixed) 


43.71/36 


17.67 


Id 


1.0 


0.98 


0.77 


9.3 X 10-« 


24.68/34 


87.92 



Table 1: The best-fit values of the parameters, sin^ 20, AF, Xb, X^^p, the Xmin' ^'^'^ ^^^ 
goodness of fit (g.o.f.) for fits to the SK1258 ES spectra. 

We find that the fits are by-and-large robust; using the separate day and night 
data or the day-night averaged data affects them only marginally. Further, the best-fit 
parameters have not changed much and they are well within the 90% CL. allowed regions 
of the previous analysis using the 1117-day data. For example, for the fit to the the day- 
night averaged data sample, in the SSM case {Xb = Xhep = 1), the best-fit values of 
sin^ 29 and AF have changed < (5 - 10)% but the quality of fit has significantly dropped 
(g.o.f. = 76.7% for the fit to the 1117-day data). When Xb and Xhep are allowed to 
vary, the best-fit value of (a) Xhep is tiny, and (b) AF is larger than earlier and closer 
to the value for the SSM fit, a welcome feature. The goodness of fit is also a bit lower. 
For the fit to the separate day and night data, the differences between the fits to the 
1258-day and 1117-day data broadly follow the same pattern as above. 



3 VEP Fits to the SNO CC and SK ES spectra 

The SNO collaboration has published the first measurement of the spectrum of electron 
kinetic energy produced from CC deuteron disintegration by solar neutrinos. The results 
are presented in eleven electron kinetic energy bins ranging from 6.5 MeV to 13.0 MeV. 
This provides a new check on solutions to the solar neutrino puzzle. The best-fit values 
of the VEP parameters obtained from fitting this data sample are shown in Table 2. 

When the SNO CC and SK 1258-day ES data (separate day and night or day-night 
averaged) are simultaneously addressed within the VEP picture, the best-fit values are 
remarkably close to each other in the SSM fits (2a), (2c), and (2e) though the goodness 
of fit is low. When Xb and Xhep are allowed to vary, then the fits to the combined SK 
ES and SNO CC spectra improve. These results are also presented in Table 2. 



Fitted 
Data 


Case 


sin^ 26 


AF 

(10-24) 


Xb 


-^hep 


xL„/d-o.f. 


g.o.f 

(%) 


SNO CC 

Spectrum 

alone 


2a 


0.70 


0.20 


1.0 
(fixed) 


1.0 
(fixed) 


9.74/9 


37.19 


2b 


0.72 


0.20 


1.09 


6.6 X 10-' 


9.70/7 


20.62 


SNO CC & SK1258 

(Day-Night averaged) 

Spectra 


2c 


0.72 


0.21 


1.0 
(fixed) 


1.0 
(fixed) 


42.04/28 


4.30 


2d 


0.66 


0.99 


0.62 


1.9 X 10-^" 


24.60/26 


54.17 


SNO CC & SK1258 

(Separate Day & Night) 

Spectra 


2e 


0.72 


0.21 


1.0 

(fixed) 


1.0 
(fixed) 


56.60/47 


15.93 


2f 


0.66 


0.99 


0.61 


9.3 X lO-^'^ 


37.79/45 


76.84 



Table 2: The best-fit values of the parameters, sm^29, AF, Xb, Xhep, the Xmin' ^'^'^ the g.o.f. 
for fits to SNO CC and SK1258 ES spectra. 

As seen in the fit to the SK ES spectra alone (Table 1), here again the fits do not 
change whether the separate day and night data or the day-night averaged data are 
used. In view of this, in the next section, we use only the day-night averaged data. 

We show in Fig. 1, the allowed regions for the parameters sin^ 2^ and AF at 90% 
C.L. obtained by fitting the SK 1258-day ES spectrum (dot-dashed line), the SNO CC 
spectrum (small-dashed line), and for a combined fit to the SNO CC and the SK ES 
spectra (solid - for the SSM case - and dotted - when Xb and Xhep are allowed to vary). 
It is noteworthy that for the SSM, the allowed regions from the different fits overlap very 
nicely. 



4 VEP Fits to total rates and the SNO CC and SK 
ES spectra 

The SK ES and SNO CC spectrum measurements are excellent filters for choosing be- 
tween alternative solutions of the solar neutrino problem. But it has to be borne in 
mind that they are both sensitive only to the highest energy *B neutrinos. Further 
constraints on any putative solution come from the total solar neutrino rates measured 
by the Chlorine and Gallium experiments. In order to be acceptable, the VEP solution 
must also be confronted with these results, to which we now turn. The data used in the 
X^ analysis of total rates are given in Table 3. 



Experiment 


Chlorine 


Gallium 


SKES 


SNO CC 


Observed Hate 
SSM Prediction 


0.337 ±0.030 


0.584 ±0.039 


0.4594 ±0.016 


0.347 ±0.029 



Table 3: The ratio of the observed solar neutrino rates to the corresponding BBP2000 SSM 
predictions used in this analysis. The results are from Refs. [^, ^, |l[.The Gallium rate is 
the weighted average of the Gallex, SAGE, and GNO results. 

The best-fit parameters obtained by fitting the total rates are presented in Table 4 
(sets 4a and 4b). Here the goodness of fits are not high. In order to ascertain how the 
total rate measurements mesh with the SK ES and SNO CC measurements, we have 
performed combined fits to the rates and spectral data. These results are also shown in 
Table 4. Note that the SSM is disfavoured by this analysis. 



Fitted 
Data 


Case 


sin^ 26 


AF 

(10-24) 


Xb 


-^hep 


xLn/d-o.f. 


g.o.f 

(%) 


Total Rates 
alone 


4a 


1.0 


1.66 


1.0 

(fixed) 


1.0 
(fixed) 


5.50/2 


6.39 


4b 


1.0 


1.66 


0.74 


1.0 

(fixed) 


2.40/1 


12.13 


SNO CC, SK1258 
ES spectra & Total 

Rates (incl. 

SNO & SK) 


4c 


0.72 


0.21 


1.0 
(fixed) 


1.0 

(fixed) 


117.62/32 


1.0 X 10-9 


4d 


1.0 


1.60 


0.77 


6.8 X 10-« 


33.54/30 


29.96 


SNO CC, SK1258 
ES spectra & Total 

Rates (excl. 

SNO & SK) 


4e 


0.72 


0.21 


1.0 
(fixed) 


1.0 
(fixed) 


92.43/30 


2.8 X 10-6 


4f 


1.0 


1.62 


0.77 


1.8 X 10-^ 


31.04/28 


31.53 



Table 4: The best-fit values of the parameters, sin^ 26*, AF, Xb, X^^p, the x^in, and the g.o.f. 
for fits to the total rates along with the SNO CC and SK ES spectra. 

When we fit the total rates and the SNO CC and SK ES spectra together, the SNO 
and SK total rates and the corresponding spectra cannot be considered to be statistically 
independent pieces of data. However, their exclusion from the fits would also not be 
entirely correct as the fitted spectral ratios (observed/SSM) do not completely determine 
the total rates. We have therefore chosen to examine two extreme alternatives; i.e., in 
the fits to the spectra and the total rates together, (a) in one case we include the SK 
and SNO total rates in the fit, and (b) in the other case we have excluded them. It 
is gratifying that the best-fit values of the parameters of the two cases (4c, 4e) or (4d, 
4f) are not very different, a reflection of the comparative energy independence of the 
(observed/SSM) spectra. 

The best-flt values of the combined total rates + spectra flts for the SSM case (4c, 
4e) are within the 90% C.L. allowed regions obtained from flts to the spectra alone (see 
Fig. 1). 
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5 SNO NC expectation, Comparison of different fits 

SNO will soon start taking data for a calorimetric measurement of neutral current (NC) 

deuteron disintegration {u + d ^ u + p + n) to which all active neutrinos contribute 

equally. A measure of neutrino oscillations is provided by the ratio of the NC and CC 

rates, Rnc and Rcc, which is somewhat less sensitive to theoretical uncertainties than 

the rates themselves. We define 

Rnc 



R 



SNO 



R. 



CC 



(6) 



For the best fit values of parameters obtained in the earlier sections, we present the 
predicted values of Rsno in Table 5. For an assessment of the various fits, for every set 
of best-fit parameters we also show in Table 5 the following quantities: the calculated 
rates for the Chlorine, Gallium, SK ES and SNO CC observations, and also the x^ 
obtained for the calculated SK ES and SNO CC electron spectra. 



Case 


CI 

(0.337 ±0.030) 


Ga 

(0.584 ±0.039) 


SK 


SNO 


Rate 
(0.4594 ±0.016) 


(spect) 


Rate 
(0.347 ±0.029) 


x' 

(spect) 


Rsno 


la 


0.516 


0.940 


0.479 


29.10 


0.353 


11.65 


1.00 


lb 


0.437 


0.796 


0.462 


11.38 


0.380 


20.96 


0.71 


Ic 


0.508 


0.940 


0.476 


29.40 


0.341 


10.73 


1.04 


Id 


0.436 


0.798 


0.458 


11.50 


0.380 


21.62 


0.71 


2a 


0.522 


0.941 


0.502 


45.50 


0.356 


9.74 


0.99 


2b 


0.523 


0.942 


0.522 


78.60 


0.356 


9.70 


1.07 


2c 


0.510 


0.940 


0.481 


29.96 


0.343 


10.35 


1.03 


2d 


0.496 


0.858 


0.461 


13.14 


0.418 


12.53 


0.51 


2e 


0.506 


0.940 


0.479 


30.25 


0.336 


10.15 


1.05 


2f 


0.495 


0.859 


0.459 


13.20 


0.416 


12.56 


0.51 


4a 


0.443 


0.609 


0.605 


355.6 


0.520 


23.18 


0.67 


4b 


0.339 


0.596 


0.448 


23.60 


0.385 


12.37 


0.67 


4c 


0.511 


0.940 


0.482 


30.26 


0.344 


10.29 


1.03 


4d 


0.352 


0.613 


0.456 


18.40 


0.400 


12.94 


0.67 


4e 


0.511 


0.940 


0.480 


29.75 


0.344 


10.44 


1.03 


4f 


0.350 


0.608 


0.456 


18.40 


0.400 


12.97 


0.67 



Table 5: The total rates for the different experiments obtained by using the best-fit values of the 
VEP parameters presented in Tables 1, 2, and 4. Also presented are the x^ for fits to the 
observed SK ES and SNO CC spectra using these values of the parameters. The predictions 
for Rsno are also shown. 

From Table 5 it can be seen that fits to the spectra by themselves give best-fit 
parameters grossly disfavoured by the total rates measured in the Chlorine and Gallium 



experiments. The most acceptable results are obtained in cases (4d, 4f) - combined 
fits to the SK ES spectra, SNO CC spectra, and the Chlorine and Gallium total rates 
measurement with Xb = 0.77 (see Table 4). These values of parameters will be further 
strengthened if Rsno is measured around 0.67. 

6 Conclusions and Discussions 

In this work, we have considered the VEP oscillation explanation of the solar neutrino 
problem within a two- flavor (active) scenario in the light of the SNO CC measurement, 
the SK 1258-day ES data, and the total rates from the Chlorine and Gallium experi- 
ments. 

We find that the VEP fits to the SK ES data are robust. The best-fit parameters 
have changed only marginally from the earlier fit to the 1117-day data. 

In the literature there are many variants of the data-fitting procedure for solar neutri- 
nos. The solution preferred by Nature should fit all the solar neutrino results collectively 
and also individually. A good fit to one (or some) of the results may do rather poorly 
when confronted with another piece of datum. A good fit to all the data taken together 
may mask bad fits to some of the results taken in isolation. This caveat is overlooked 
in much of the current literature. We have attempted to fit the SK ES and SNO CC 
spectra individually and jointly, and also the total rates by themselves and jointly with 
the spectral data in an attempt to find a region of parameter space that is acceptable 
in every respect. These results are given in Table 5 and in Tables 1, 2, and 4. We find 
that the best-fit parameters of the cases (4d) and (4f) (see Tables 4 & 5) best meet 
the requirements alluded to above. Note that these fits correspond to Xb ~ 0.77 and 
predict Rsno ~ 0.67. 

During the passage of the neutrinos from their point of production to the solar 
surface, interactions with the ambient matter, responsible for the MSW effect, become 
important. In the presence of VEP and this MSW contribution, the effective neutrino 
mass matrix in flavor space takes the form 



M = i 

2 



E^AF cos 29 - 2V2GFne{r) E^AF sin 29 
E^AF sin 29 -E^AF cos 29 



(7) 



where we have dropped an irrelevant part proportional to the identity matrix. The 
MSW contribution in (|^ inside the sun turns out to be several orders of magnitude 
larger than the terms due to VEP that we have discussed in this work. Recall that we 
have assumed the neutrinos to be massless. For maximal mixing (e.g., the preferred 4d 
and 4f cases), there is no resonance effect and, in fact, till such time that the neutrino 
emerges from the sun, the MSW contribution controls the masses in (0). Inside the sun, 
the Ve is, therefore, a mass eigenstate to a very good approximation. The effect of VEP 
oscillations begins to manifest itself only from then onwards. 



The VEP oscillation wavelengths favored by the data are comparable to the distance 
of the sun to the earth. Therefore, seasonal effects are expected in this picture. When 
more seasonal data from SK accumulates, it will provide a good check on this scenario. 

Since the publication of the first results by the SNO collaboration, a number of 
analyses of the data in the mass-mixing vacuum oscillation and MSW approaches have 
appeared p. The quality of these fits are comparable to those found in the VEP 
picture. Here too, the parameters which best fit the rates alone are somewhat different 
from those that best fit the rates together with the spectra pf. Further sharpening of 
the data from new and ongoing experiments, it is hoped, will help to distinguish between 
these alternatives. 
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Figure 1: The 90% C.L. allowed regions in the sin 29 - AF plane. The area enclosed by the 
solid (dotted) line is allowed by the SNO CC + SK1258 ES spectral data for the SSM (the 
model where Xb and Xh^p are allowed to vary). The area allowed by the SNO CC (SK1258 
ES) spectrum alone for the SSM is enclosed by the small dashed (dot-dashed) line. The best 
fit points have been indicated. 



11 



